Assisted tunneling of a wave packet between square barriers 
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The assisted tunneling of a wave packet between square one dimensional barriers is treated ana- 
lytically. The survival probability is calculated exactly for a potential mimicking a constant electric 
field with arbitrary time dependence. The pole spectrum of the properly normalized unperturbed 
wavefunctions determines the decay time. The tunneling probability is enhanced by the perturba- 
tion. The results are exemplified for a simplified model of nuclear a decay. 
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I. INTRODUCTION 



The quantum tunneling paradigm explains diverse pro- 
cesses, as currents in Josephson junctions and nuclear al- 
pha decay. Energetically forbidden regions for a classical 
object, become accessible for a quantum system. 

Gurney and Condon and Gamow[3, 0] , used the tun- 
neling model to reproduce the experimentally measured 
huge range of a decay lifetimes of similar unstable nu- 
clei. The scheme has withstood continuous scrutiny in 
the intervening years [1, 0] • 

In a previous work we investigated the tunneling of a 
one dimensional metastable state between delta function 
barriers, excited by a time dependent potential. [l| By in- 
specting the modification of the spectrum of the poles 
ruling the decay process, we estimated a non-negligible 
enhancement of the tunneling probability when an exter- 
nal time dependent perturbation acts on the system. 

In the present work we find exact solutions to the as- 
sisted tunneling process, with square barrier potentials 
instead of delta functions. Square potentials provide a 
more realistic description of the barriers involved in de- 
cay processes. 

Decay times are enormously different form the nat- 
ural time scales of nuclear and even atomic phenom- 
ena. It is not possible to follow the decay process nu- 
merically. Analytical expressions are then of the utmost 
importance 11, ^, £]. The method and formulas presented 
here aim at contributing in this direction. In the next 
section we present the model and the analytical solution 
method. In section 3 we solve the equations for a time 
harmonic, linear potential and evaluate the effect of the 
perturbation on the tunneling process. Some mathemat- 
ical details are presented in appendix A. 



II. ANALYTICAL METHOD FOR ASSISTED 
TUNNELING 

The Schrodinger equation for a one dimensional system 
between square barriers is 
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dt 2 TO dx^ 



+ j {Q{d^\x\)e{\x\-xo))^ (1) 



In the following all our units are either fm for length 
and time, and fm"^ for energies, momenta and potential 
strength 7. 

We consider an initial Gaussian wavepacket located t — 
in the region between the barriers. 



= 0) = e" 
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where A is the width parameter of the packet. This 
packet is a reasonable model for the metastable state and 
facilitates the analytical evaluation of the decay ampli- 
tudes. However, the method presented here is not limited 
to gaussian packets. 

Inside the barrier region the packet disperses if 
A < xq. , but, soon enough swelling comes to a halt due 
to the presence of the barriers. It can then spread only 
through the tunneling process governed by the barriers. 

The even Xeix) — —j^ — jj-r and odd Xo{x') 

ipo{k) 



/7rne(fc) 

stationary states of eq.([T]) for energies below 



the barrier height 7 are 
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where k is defined in eq.([7]), and sign{x) 
x> Q,x <Q. 
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where we have extracted a factor of \pK from the normal- 
izations for convenience. 

The set of even-odd functions is orthonormal and 
complete. ^Sl, The zeros of the normalization factors 
govern the exponential decay of the metastable wave 
functions in the inner region[l|. 

The normalization factors are given by 
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63 = 6" (^»+'') 
61 = 6^ ^° 

S2 — sin(k d) 
Si = sin{k xq) 
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We assume a time dependent spatially linear pertur- 
bation potential 



V{x,t) =nx G[t) 



(8) 



with /X a coupling constant. For a spatially constant time- 
harmonic electric field of intensity Eq interacting with an 
a particle of charge 2 \e\, we have 



H G{t) = 2 |6| Ea sin{ujt) 



(9) 



In order to eliminate the explicit space dependence of 
the perturbation term in the Schrodinger equation, we 
apply the unitary transformation 



a = X Q -\- I — at 
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G(t) dt 



(10) 



The Schrodinger equation([T]) including the perturbation 
of eq.® becomes 
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The wave function $ is expanded in the complete set 
of even and odd states of the unperturbed Schrodinger 
equation ([T]) including states with energies above the bar- 
rier 
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The amplitudes ee,o consequently obey the time evolu- 
tion 
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where we have used the superposition integral calculated 
in appendix A 
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and a dot denotes a time derivative. 

The exact solution to eqs. (IT5)l for an initially symmet- 
ric wave packet such as the one of eq.© is 
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where Ce{k, 0) corresponds a well localized A << xq wave 
packet of the form of eq. ^ . Other initial wave functions 
merely change the functional dependence of Ce(fc, 0). 



III. ASSISTED TUNNELING SOLUTION 

We are interested in the decay of the wave packet from 
the internal zone leaking out to infinity. As discussed 
in the most important contribution to the integral 
in eq. p^ originates from the poles in the amplitudes of 
eq. (lisp . This is indeed correct for the unperturbed decay 
process. However, for the perturbed process, the ampli- 
tudes are oscillatory. The poles lie in extremely close 
proximity to the real momentum axis. When the inte- 
gration over this variable is effected, the argument of the 
harmonic functions in cq. (jl5p can oscillate wildly, sup- 
pressing the contribution of the poles, (see figures 2 and 
3 below) However, if the argument of the harmonic func- 
tions is of order one, the poles will still dominate. This 
can be achieved by choosing a sufficiently high frequency 
for the external perturbation of eq.®. 

The normalization factors around their minima can be 
cast in the form 



Aj,(e,o) {k'' - fcj,(e,o)^)^ + f3j,{e,o) (16) 



where j enumerates the zeros of n,,, Uq 
The minima are located close to 



{2n+l) 77 p 

2 {l+pxo) 
n t: p 
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p = y^2~mj 



(17) 



Eas. (|17p provide starting values for the numerical search 
of the location of the poles of the wave function. The 
poles in the complex momentum plane, appear in pairs 
located symmetrically above and below the real momen- 
tum axis. The stationary wave functions of eqs. (|3l4p con- 
sist of incoming and outgoing parts. The standard S 
matrix, reflecting the behavior of the outgoing boundary 
condition at spatial infinity, has poles in the lower half of 
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FIG. 1: 10"^°ne(fc)~\ and 10"^°no(fc)"^ , as a function of k 
in units of fm~^ for a typical a decay barrier 



the complex momentum plane, whereas its inverse, cor- 
responding to the incoming condition, has poles in the 
upper half. Consequently, both sets of poles in the up- 
per and lower complex momentum plane appear when 
incoming and outgoing pieces are present. TC|. 

As the poles are extremely sharp in position, it was 
necessary to use a quadruple precision algorithm on an a 
cluster mainframe to find them. The imaginary parts of 
the poles are typically many orders of magnitude smaller 
than the real parts. 

Figure 1 shows 10"^"ne(fc)"\ and 10^^"no{ky^ , for 
the parameters m = 3727Mey — 18.88/to~^, = 
12fm,d= 22/m,7 = 22MeV = 0.11/to"\ correspond- 
ing to a typical a decay barrier. The spikes are due to 
the extreme closeness of the minima of the normalization 
factors to their complex zeros. 

When eq. (fT2|) is evaluated by contour integration in 
the complex k^ plane, the contour has to be closed from 

below. In the lower half-plane the convergence is insured 

_ ■ fc^ t 

by the exponential e ' ^ ^ . xhe poles are separated from 
each other and the integrand drops essentially to zero 
outside the pole, hence, the wave function consists of 
an incoherent sum of the residues of the different poles. 
The negative imaginary part of each pole induces a time 
decaying exponential. Each exponent d_etermines a dif- 
ferent decay constant and decay time.|l,] If the original 
wave function is even in space, as in eq.©, only even 
poles show up in the unperturbed wave function. 
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The argument of the harmonic functions in eq. (jlSp in- 
cludes a potentially large denominator near the minima 
of the normalization factors. The product of these even- 
odd factors near a minimum is depicted in figures 2 and 
3. This product is not only finite, but of order one. We 
were not able to produce an analytical formula for this 
product, although the numerical values are very sugges- 
tive. 
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FIG. 2: no{kj,e) nj^e ^ , around the minima of ne{k). Param- 
eters as in figure 1 
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FIG. 3: no{kj,e) rij^e ^ , around the minima of no{k). Param- 
eters as in figure 1 
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(19) 



Figure 2 depicts {ne{k) rio(fc))^^, around the minima 
of ne{k). Note the degree of accuracy needed for the 
location of the minima. The parameters are those of 
figure 1. Figure 3 depicts (ne(fc) no{k))^^ around the 
minima of no{k). 

Expanding the harmonic functions in cq. (jl5p , the exact 
expression for the wave function in the region between the 
barriers becomes 



^{x,t) = <^^e{x,t) +^oo{x,t) +^eo{x,t) +^oe{x,t) 

(18) 



where 
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The series in eas. (|19l20l21l22p can be expressed in terms 
of standard functions 



Jo(uj>)Hi(uj^o) - Ji('Uj,o)Ho(wj,o 



Jo(%\o)Hi(lij,o) - Ji(Wj,o)Ho(Wj,o) 



2 Jo(Mj,e) 
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Jo(Uj,e)Hi(?ij,e) - Ji(Uj^e)Ho(Uj^e) 



(23) 



where Jo,i denotes the Bessel function of the order (0,1), 
and Ho.i is the Struve H function. 11] 

For long times we can use the asymptotic expansions 
of the Bessel and Struve functions to obtain llli 



^ee{x,t) ^ Gj{x,t) 



^'eo(a;,t) 
*oe(a;,t) 



J2 H,{x,t) f^ + Oiuj} 



+ 0{u^J)] (24) 



From eq. (p4)) it is clear that for long times the even-even 
component is essentially the dominant one. We can now 
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FIG. 4: S(t) of ea. (|25p . without, full line, and with harmonic 
perturbation, dashed line, see text for parameters 

evaluate the influence of the perturbation on the decay 
time. Figure 4 depicts the nonescape probability of the 
system in the region between the barriers 



S{t) = 



\-^{x,t)\^ dx 



(25) 



for the square barrier case, we have used the parameters 
relevant for a decay mentioned below eq. ([T]) and an intial 
wave packet of width A = 4/to. The lower curve corre- 
sponds to (see eq.®) an electric field of = O.lVolt/m, 
and a frequency of w = IOtt lO^"^ Hz. The aim of this work 
is to present the analytical solutions and show the accel- 
eration of the decay process. A more realistic choice of 
parameters demands a thorough numerical evaluation of 
the integrals outside the poles and is postponed for future 
work. For a time harmonic perturbation®, there are two 
major corrections to the unperturbed decay rate. The 
amplitude of the wave function is to lowest order multi- 

plied by a factor of the form (1 — ct^/4), with c ~ — ^ — ^ 



from the series in 



This factor depletes the wave at 



times of the order of the period of the harmonic pertur- 
bation for moderate values of /i, much faster than the 
actual decay time determined by the poles. A second 
correction already noted in[l|, arises when considering 
the higher order poles j > 1 in the sum of eq. HH) . The 
latter generates a different decay time for each pole[l[. 
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Figure 5 shows the results for delta function barriers 
V 8{x—X{^-\-vb(x-\-XQ)^ with parameters corresponding to 
the square barrier case v = ^{d — x^)- The unperturbed 
decay time of around lO^^^sec is orders of magnitude 
smaller than the square barrier one. The assisted tunnel- 
ing is barely noticeable. 

In summary, we presented a simple analytical method 
to calculate assisted tunneling. The technique can be 
generalized to three dimensions, and to other potentials 
and initial metastable states. Such endeavor is worth 
pursuing due to the dramatic impact assisted tunneling 
may have on nuclear technologies as well as on solid state 
devices. The question of electron screening for the nu- 
clear decay acceleration must also be addressed to pro- 
duce experimentally meaningful predictions. 
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IV. APPENDIX A: PERTURBATION MATRIX 
ELEMENT 

The matrix element of eq. does not vanish for tran- 
sitions between even and odd states solely. The matrix 
element evidently diverges and we seek to evaluate the 
divergent part. The most important contribution comes 
from the outer region of jxj > d. In the inner region there 



will be a finite contribution that is negligible as compared 
to the S function contribution of the outer region. Conse- 
quently we can use the outer region expression through- 
out. 

The functions Xe,o are then replaced for all x by 
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(26) 



Taking advantadge of the symmetry properties of the 
wave functions, the matrix element of ea. (|14p becomes 
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^ , , . , / ^, ''^'^ ■'^ A lengthy but straightforward calculation then gives the 

i{-k + k') result of eq.IlM]) 

-i(k+k')x 

M.) = -^'^*i^*Z-(^ (27) 

where the tilde denotes a factor that depends on Jc' in- k , 

stead of k. The all important contribution comes from -^fc',fc = ^ j-^-j ^ (k)^^^ ~ ^ ^ ^'^^^ 

k' ^ k 

M,,.j, « 2 ^ - j^*) '^Ifr 11 (28) 

ne(k)no[k') 



